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Abstract 

We present evidence in the Schwinger model that rooted staggered fermions may correctly 
describe the m < sector of a theory with an odd number of flavors. We point out that in 
- --_ QCD-type theories with a complex- valued quark mass every non-chiral action essentially 

^-^ • "borrows" knowledge about the ^-transformation properties from the overlap action. 

o' 

The staggered fermion operator D^^ describes Nt = 2^^/^ flavors ("tastes") of quarks in d 
dimensions. In order to obtain a single quark flavor in dynamical fermion simulations, it has 
become customary to include the A^^-th root of the staggered determinant as a weight factor. 

In a recent preprint Q Creutz has argued that this procedure makes it difficult to obtain the 

right continuum limit. His argument is based on the observation that the staggered determinant 

^ ' is strictly a function of m^, whereas the determinant of a true single flavor theory also contains 

•^ odd powers of m, arising from the unpaired chiral modes of the Dirac operator. 

Q In a reply P] Bernard et al. have noted that the rooted determinant always corresponds 

■^ , to a theory with positive fermion mass, regardless of the sign of the staggered mass term. In 
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addition, they point out that rooting is a non-analytic procedure. Accordingly, they discuss 



O ' how taking the continuum limit flrst can lead to odd powers of |r7i|. 

In this note we want to address how two observables, the chiral condensate and the topo- 
logical susceptibility, behave for non-positive quark masses. We demonstrate in the Schwinger 
Qh! model [3, that even with rooted staggered fermions these quantities can be measured for m < 0, 
pH ' by simply reinserting the phase information which gets lost by starting with a doubled action 
IJ and taking the square-root. This "recipe" can be justifled by reference to the overlap action, 
but it is beyond the realm of a purely Lagrangian approach with staggered fermions. 

The numerical data presented below are based on a sample of 1000 gauge conflgurations at 
F^ , P = 9.8 on a A^xA^ = 28^ lattice. We compare (rooted) staggered to overlap fermions [4J, both 
actions deflned with one step of APE smearing. For details see 1310]. Throughout, we use the 
coupling e to deflne the physical scale and set the lattice spacing to a = 1. 

First, we consider the scalar condensate. We are interested in its continuum limit in a flxed 
physical volume (cf. |6;), and for Nf = 1 its presence is due to the anomaly and does thus 
not signal spontaneous symmetry breaking. For overlap fermions and m > its bare version is 

Xscai_v^ (det(D-)^/E^^) - 1 

e m (det(D-)^/) ' i-^ ^^ 

where the sum runs over all eigenvalues of the massless operator 1)°^, except for the doublers (on 
topologically non-trivial backgrounds) at 2p. Here, p is a parameter in the overlap construction 
which in our runs is set to 1. The measure in (Q) is quenched with an explicit weight factor 

rn 

det{D'^) = Yl{il--)X + m). (2) 



For staggered fermions and m > the definition of tlie bare condensate is analogous 
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XZ^ ^ v/? (det(Dy//^E^) 
e 2m (det(L'^)^//2) ^ ^ 

except tliat tlie sum over all eigenvalues plus the prefactor 2~^ = Nt^'^ provide a taste projection 
in the valence sector and the ensemble average is taken with the square-root of the determinant 

detiD:^)=l[i\ + m) (4) 

to correct for the extra degree of freedom in the sea sector. For m > the condensates (pQ Ej) 
have a logarithmically divergent contribution, and after subtracting the latter, they seem to 
yield a universal quantity in the continuum, regardless whether Nf = 1 or Nf = 2 (where there 
is no rooting issue in 2D) jHj. Below, this subtraction is always implicit. 

We now evaluate ((H ^ for arbitrary real m with details as specified in t^. Triple curves 
in Fig.^ indicate la error bands. Errors are correlated, since the same gauge configurations 
have been used. We first discuss the Nj = 2 case. Here, our overlap and staggered data are 
in good agreement. They are both odd in m, as expected with an even partition function and 
one derivative w.r.t. the quark mass. The Nf = 1 case is more subtle. Here, our overlap data 
are smooth near m = 0, reproducing the analytic result Xscai/e = exp(7)/(27r^/^) = 0.1599... 
by Schwinger j3]. On the other hand, the rooted staggered data pass through the origin, as 
they have to, since their partition function is even in m, and a derivative w.r.t. the quark mass 
yields an odd function. The 1-flavor overlap curve is neither even nor odd, since the operator 
being odd is not sufficient for inducing this property in the observable. The functional measure 
is even in m with a straightforward ("naive") evaluation of (JHIEI) in case of staggered fermions. 
It is, however, not even with an odd number of overlap fermions, as is obvious from the fact 
that on a background with index u the overlap determinant has the same sign as m^-fl'^L 

With this information at hand, it is now easy to understand how the failure of the "naive" 
staggered rooting procedure at m < could be corrected. What is missing is the same informa- 
tion that gets lost if one first squares the overlap determinant and then takes the square-root. 
This suggests that - rather than naively evaluating the product in Q for any m - one should 
define the functional determinant of a rooted staggered fermion at negative mass via 

det{D^\Nj = l,-|m|) ^ (-l)''^' ciet(Dfl[)^/^* (5) 

in close analogy to det{D'^,^ = (— l)'*^' det(DP^|) (see below). We have tried the definition 
(0), together with a sign-fiip in the condensate (see below) for A^ = 1; the result is shown with 
dash-dotted lines in Fig.[TJ This definition lets the "smart" staggered action reproduce the 
overlap answer for the 1-flavor condensate at (sufficiently) negative mass. Together with the 
standard definition at positive mass, the staggered answer is thus qualitatively correct, except 
for a small dip at m = 0. This dip gets narrowed (horizontally), as the coupling tends to zero 
(/3 — >■ cx)), with details given in [6^. Still, since the value at m = remains an exact zero at any 
j3, the continuum extrapolated staggered curve is discontinuous at this point. This shows that 
in general the limits m ^ and a ^ do not commute with staggered fermions [HI El- 

Next, we consider the topological susceptibility in the same finite box. It is defined through 

Xtop _ P {det{Dm)^f z/2) 



e2 m {det{Dm)^f) 
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Figure 1: Scalar condensate in the A^ = 1, 2 dynamical theory versus quark mass. For Nf = 1 
and 171 = the overlap condensate successfully reproduces the analytic value by Schwinger. For 
Nf = 1 and m < only the "smart" staggered curve [with an explicit factor (—1)''] is correct. 
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Figure 2: Topological susceptibility in the A^ = 1, 2 dynamical theory versus quark mass. For 
Nf = 1 and m < only the "smart" staggered curve [with an explicit factor (—1)'^] is correct. 
The quenched topological susceptibility is included to indicate the asymptotic m -^ +cxd value. 



and we choose to evaluate the index v with the overlap operator, but other options would 
just bring different cut-off effects. The determinant in (jHI) is meant to represent the eigenvalue 
product on the r.h.s. of Q for overlap fermions. For staggered quarks there is the "naive" 
option to use the r.h.s. of ^ for all m, versus the "smart" one to use it only for m > 0, and to 
continue to m < via (0). 

Our results for the topological susceptibility with real m are shown in Fig.|21 For Nf = 2 
both actions are in perfect agreement; in particular both curves are strictly even in m. For 
Nf = 1 the overlap curve goes right through the origin, but it is not exactly odd in m. It is 
negative for m < 0, in spite of the observable u^ being positive or zero, simply because the 
determinant may have either sign. By contrast, the "naive" staggered curve is still exactly even 
and positive. Being smooth, it must show a residual susceptibility at m = 0, a pure cut-off 
effect which tends to zero, under (3 — > cx), with a rather high power of ae 0. Again, defining 
a "smart" staggered action at m < via (0) basically flips the sign of the staggered 1-flavor 
curve (for small enough \m\) and makes it agree with the overlap data. Still, the derivative 
from the right, evaluated at m = 0, is a "malicious" observable where the staggered continuum 
limit is not correct. It is always zero, whereas the overlap counterpart indicates that the true 
continuum- limit is non-vanishing. This is another example of the non-commutativity of the 
771 —i> and a — *> limits in the Nf = 1 theory jB] . Regarding the m > piece of the curve, we 
find it most amusing that the staggered version consists of a rather large tower of even (in m) 
contributions which - all together - manage to mimic an almost linear behavior. 

We add a remark regarding m < 0. Our graphs in Figs.^Elare cut off at m/e = —0.1. This 
is because the sign problem in the 1-flavor theory hits us hard to the left of this point, and the 
error-bars explode. Note that there is no sign-problem in the 2-fiavor theory, regardless of the 
action, and the "naive" staggered 1-fiavor determinant does thus not induce a sign problem. 
By contrast, the "smart" staggered 1-fiavor theory has the same sign problem as the overlap 
version, and this speaks in favor of (0) being an honest approach to the m < region. 

We now extend the discussion to QCD with arbitrary complex quark mass, assuming that 
the fundamental ^-term in the action is zero (if not, it would entail a trivial generalization). 

In the continuum, a complex quark mass m = \m\e^^ can be rotated into such a 6'-term, and 
the derivation builds on chiral symmetry jHj. In a lattice regulated theory the same holds true 
if the lattice action respects exact chiral symmetry. To this end it is useful to note that the 
standard definition of the massive overlap operator, D'^ = (1 — y)D°^ + m, is not adequate 
for complex m. For reasons to become obvious, we represent the action in the general case as 

SZ = ^+^°>+ + ^-/^°>- + m^_^+ + m*ilj+ip- (7) 

(compare jH] and P) where the positive/negative chirality fields are defined through 

^± = P±7/>, iP± = ijP^, (8) 

^± = ^(1±75), P± = ^(l±75), 75 = 75(21-1), i = l-^D°- (9) 

and hat/tilde-quantities implicitly depend on the gauge field. By defining rotated fields 

V^; = e±'^/2^±, ^P'^ = ^P^e^''/' (10) 

the overlap action ((7j) can be brought into the form 

{s^y = ^;/^°>; + ^^/^°>^ + |m|^^^; + \m\tp'_,^'_ (ii) 
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which reflects only the modulus of the quark mass. In addition, one has to take into account 
the anomalous transformation of the 1-flavor path-integral measure [TC 



DipDiP — > Dip' Dip' = Dip Dip detf f- ,^ ) (12) 

^ dip dip '^ 

and with Sip' /Sip = e^"'^^^'^, Sip' /Sip = e^^^^l"^ and the index theorem JT] it follows that 
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Sip Sip ' 
With a chirally symmetric action one can thus rotate the effect of a complex mass into a term 

Se = -ive (14) 

in the action and stay with a positive fermion mass \m\. In a theory with Nf degenerate 
overlap fermions it is convenient to parameterize each mass as m = \m\e^^^^f ^ since this choice 
again yields the 6'-term fll4j) . It is straightforward to generalize the discussion to an arbitrary 
complex Nf X Nf mass matrix, and the result is the same if one lets 6 represent the phase of 
the determinant of this matrix i^Sj. The important consequence of the ^-term (I14j) reflecting the 
anomalous transformation behavior of the QCD path-integral measure is that the full theory 
is periodic under a simultaneous 2tt /Nf shift in the phase of each quark mass. 

With staggered fermions it makes no sense to "naively" evaluate the action with a complex 
mass. However, the discussion in the overlap case has shown that it is sufficient to have an 
explicit 6'-term (I14j) and a fermion action with a positive mass. Hence, if we assume that the 
Nf/Nt-th root of the staggered determinant yields the right continuum limit for Nf positive- 
mass fermions, then it is perfectly reasonable to stay with that root and to introduce an explicit 
^-term fll4|) in the action. The pertinent 2tc /Nf period under a common shift in the phase of 
each quark mass in the resulting continuum theory is no further assumption. 

In order to compute observables for arbitrary complex mass via the 6'-term prescription, one 
has to perform a chiral rotation in the valence sector, too. We first discuss the overlap case. 
Representing the mass term in (|7j) as \m\ip{Pj^e^^ Pj^ + P^e^^^ PS)ip and using the identities 

P+P+ = ^(l+75)(l-^/^°^) = ^(1+75)1 (15) 

P-P- = ^(l-75)(l-^^°^) = ^(1-75)1 (16) 

the 1-fiavor action with a complex mass m = \m\e^^ can be written in the compact form 

S:: = iPD°^iP, D- ^ D- + |m|e'^^^i . (17) 

In the theory with complex mass m = \m\e^^ and vanishing 6'-term, the (unsubtracted and 
unrenormalized) scalar /pseudoscalar condensate is given by 

Xscal = -{iPilP) = -{lP{P+P+ + P-P-)lP) = -{lPj++lP+4^-) (18) 

Xpscu = -{lPlillP) = -{lP{P+P+-P-P-)lP) = -{lPj+-lP+^-) (19) 

if one includes the determinant of the massive operator (fTTj) into the functional measure. We 
now express them in terms of the primed fields which refer to the theory with positive mass 

Xscai = -(V^-V'+ + ^+V^-) = -{'ip'-e~'^ip'+ + ^'+e'^'__) = cos(^) xLi - isin(^) Xpseu 

Xpseu = -{^-■ip+ - ^+^-) = -W-e~'^ip'+ - ^P'+e'^^'_) = cos(^) Xpscu - isin(^) xLi 

5 



where the latter quantities adopt a simple form which shows invariance under chiral rotations 

xLi = mi{D- + \m\l)-']) = ^(E' J^) . xUu = (Tr[75ip°^ + |m|i)-i]) = ^^ • (20) 

In these formulas the functional measure is still the one in the theory with the complex mass. 
Our previous exercise tells us that this is equivalent to the quenched measure with det(L'P^|) 
and an additional factor e'^'^^. In other words, we need to calculate the expressions in (|2(jp times 
some reweighting factors on quenched backgrounds. In a theory with Nf degenerate flavors of 
mass m = \m\e^^^^f the scalar/pseudoscalar condensate (per flavor) is thus given by 

o. ^ 1 (det(D|--|)^/ e-^ [cos(g/iV,) ^ x^ - isin(g/iV,)^]) 
^scai y (det(Z}|^^|)^/ e"^^) ^ ' 

.. ^ 1 (det(D|--|)^/ e-^ [cos(g/iV,)^ - i sin(g/iV,) E' x^) 

where the sum is over all modes of the massless theory, except for the chiral doublers, and A 
relates to A as indicated in (Q). With staggered fermions there is no alternative to defining the 
scalar/pseudoscalar condensate (per taste) through the expression 

.t ^ 1 (det(D|-^|)^//^- e'-^ [cos{e/Nf)Nr'Ej^^-iMd/Nf)j^]) 
^scai y (det(Df;;,|)^//^' e'-o) ^ ' 

with the sum running over all modes of the massless theory and the factor Nt~^ taking care 
of the degeneracy. The case of the topological susceptibility is simpler. Since the topological 
charge is invariant under a chiral rotation, we have in the overlap case the identity 



1 (det(D|^^|)^/ e'^' z/2) 
Vld^tCDJ^F^ 



^top T/ /J„j-/ nov NAT, „\vQ\ \^> 



and the analogous definition for staggered fermions reads 

1 (det(Df* i)^//^' e'"^ J/2V 
^*°P V (det(Df;;|)^//^' e-«) ■ ^ ' 

In all these expressions one may replace e"^^ -^ cos(z/6'), isin(i/^), depending on whether the 
remainder is even or odd in v (due to CP symmetry in the Yang-Mills theory). We stress 
that an attempt to directly implement (jTj) with Wilson fermions is conceptually not very much 
different from the approach where one "borrows" expressions (|22l 0^ from an action with exact 
flavor and chiral symmetry, as we have advocated for staggered fermions. 

In order to test how this reduction works in practice we return to the Nf = 1 Schwinger 
model. We restrict ourselves to 6 = ±7r/2, i.e. to the imaginary mass axis. Here, the scalar con- 
densate is real, and the pseudoscalar one is purely imaginary; Xscab Xpseu/i are shown in Fig. El 
For Xscab the overlap expression (j22p shows a smooth behavior, reproducing the Schwinger 
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Figure 3: Scalar condensate and imaginary part of the pseudoscalar condensate versus imagi- 
nary mass in the Nf = 1 dynamical theory, with the definition ()23| 0^ in the staggered case. 
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Figure 4: Topological susceptibility versus imaginary mass in the Nf = 1 dynamical theory. In 
the staggered case the definition (j^ has been used. 



value in the m — >■ limit. The staggered prescription ()24|) seems consistent, except for a region 
around m = 0, where the non-commutativity phenomenon again occurs. The absence of an 
exact matching of the zero-modes in the determinant and in the observable leads to a blow-up 
in the staggered answer, if the chiral limit is reached from the imaginary-mass direction. In 
Xpseu/i, the deviation between the two actions seems to be a mild 0{a^) cut-off effect. 

Our data for the topological susceptibility in the A^ = 1 Schwinger model with quark 
mass m = |m|e'^ and 6 = ±7r/2 are shown in Fig.|3 In this case the situation looks entirely 
inconspicuous, and we expect that there is a common continuum limit. 

We have demonstrated in a simple toy model that one may obtain sensible physics at 
negative quark mass for an odd number of staggered continuum flavors, if one uses a "smart" 
definition of the rooted staggered determinant which takes into account the phase factor that 
got lost in the rooting procedure. Without the advent of exact lattice chiral symmetry, this 
would be hard to justify, since with staggered fermions it is definitely beyond what a purely 
Lagrangian approach can do. We close with three conjectures on staggered fermions: 

1. For m > 0, taste-projected correlators built from D~^^^ and det(Dst,m)^^''^* in concert 
yield the correct continuum limit for rf = 2,4, arbitrary gauge group and Nf. 

2. For m = 0, taste-projected correlators built from D~^\ and det(Dst,o)^^''^* define a theory 
in the wrong universality class (i.e. at least one observable assumes an incorrect continuum 
value) for d = 2,4, arbitrary gauge group and Nf. 

3. The "smart" staggered determinant which effectively shifts the ^-dependence of the theory 
with m = \m\e'^^^^f into an explicit ^-term fll4j) yields the correct 6'-dependence if and 
only if the staggered continuum limit with real positive mass |r7i| is correct. 
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